Bound lattice complexes through distance selective two-body loss 
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We show that in a gas of ultra cold atoms distance selective two-body loss can be engineered via 
the resonant laser excitation of atom pairs to interacting electronic states. In an optical lattice this 
leads to a dissipative Master equation dynamics with Lindblad jump operators that annihilate atom 
pairs with a specific interparticle distance. In conjunction with coherent hopping between lattice 
sites this unusual dissipation mechanism leads to the formation of coherent long-lived complexes 
that can even exhibit an internal level structure which is strongly coupled to their external motion. 
We analyze this counterintuitive phenomenon in detail in a system of hard-core bosons. While 
current research has established that dissipation in general can lead to the emergence of coherent 
features in many-body systems our work shows that strong non-local dissipation can effectuate a 
binding mechanism for particles. 

PACS numbers: 37.10. Jk, 32.80. Ee, 34.20.Cf 



Dissipation and loss are intuitively considered harm- 
ful for any phenomenon relying on quantum coherence. 
Contrary to this it has been recognized recently that dis- 
sipation can in fact foster a host of coherent phenomena. 
Lattice gases of ultracold atoms [T] provide an experi- 
mental platform where these effects can be analyzed and 
studied in great detail. It has been theoretically shown 
that these setups permit the engineering of dissipation 
through a careful design of the coupling of the system 
degrees of freedom to a fictitious heat bath [31 13] . In this 
context it was demonstrated that a specifically tailored 
dissipation can give rise to the formation of an atomic 
superfluid [4], or BCS-like fermion pair states with d- 
wave symmetry [S] - paradigmatic examples for a coher- 
ent state of matter. Beyond this engineered case also 
dissipation stemming from local particle loss through in- 
elastic collisions can induce correlations in quantum gases 
[H [7] . Experimentally this has been shown by Syassen 
et al. for a lattice gas of molecules with strong local 
(on-site) two-body loss [S] . Here the phenomenon under- 
lying the formation of correlations is the quantum Zeno 
effect which prevents the simultaneous occupation of lat- 
tice sites by two bosons 9J and thus gives rise to an ef- 
fective hard core interaction. 

In this work we are interested in a scenario in which 
particle loss is taking place non-locally, i.e., we investi- 
gate a situation where atoms are expelled from the sys- 
tem only when occupying distant lattice sites. We show 
that in the limit of strong loss such distance selective dis- 
sipation leads to the formation of various types of long- 
lived lattice complexes that are formed by two or more 
atoms. Our work not only demonstrates that strong non- 
local dissipation features an effective binding mechanism 
for particles, but, in addition, shows that these dissi- 
patively bound complexes may even exhibit an internal 
level structure, which is strongly coupled to their external 
motion. We illustrate how this particular kind of dissi- 
pation can be achieved in practice by the resonant laser 
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FIG. 1. (a) Level scheme of pair of atoms laser-excited to 
Rydberg states. The pair state |ee) is excited via an effective 
two-photon resonance from the state \gg) when the distance 
between the two atoms is equal to the critical distance r^. Me- 
chanical forces due to the strongly repulse Rydberg-Rydberg 
interaction and radiative decay (at total rate 2rryd for two 
excited atoms) lead to loss of the atom pair, (b) In an optical 
lattice this leads to a distance selective loss of atoms. Shown 
is a pair of atoms that is created at the critical distance after 
a hopping event, leading to a two-body loss. 



excitation of strongly interacting atomic Rydberg pair 
states. On the one hand, our work thus generalizes cur- 
rent ideas related to the creation of effective long-range 
interactions induced by the laser excitation of atomic Ry- 
dberg states [I0HI3] and contributes, on the other hand, 
to recent efforts of using the exaggerated properties of 
Rydberg states to tailor and control dissipation in ultra- 
cold gases [m [15] . 

The setup we are considering consists of a one- 
dimensional lattice gas of ultracold bosonic atoms. We 
assume a deep optical lattice potential so that the mo- 
tion of the atoms is confined to the lowest Bloch band. 
The atoms can move through the lattice by tunneling 
between neighboring lattice sites with a rate J and ex- 
perience an on-site interaction U , whenever two of them 
occupy the same site. For the purpose of this work we 
assume that the on-site interaction is large, {U » J) so 
that the coherent motional dynamics of the system is well 
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captured by that of a Tonks-Giradeau gas [T5] with the 
Hamiltonian H = JYl,j=iW~j'^J+i + "'/"'j+i] US- Here, 

denotes the total number of lattice sites, and cr^ are 
the conventional (spin-1/2) lowering and raising opera- 
tors related to site j, where an up (down) spin represents 
an (un)occupied site. 

Let us now briefly outline how distance-selective loss 
can indeed be generated - more details will be given to- 
wards the end of the manuscript. The idea is to ap- 
ply a laser field with Rabi frequency 17 and detuning A 
that couples the electronic ground state \g) of an atom 
to an highly excited (Rydberg) state |e). When two 
atoms are excited simultaneously to the Rydberg state 
(i.e. they are in the pair state |ee)) they strongly inter- 
act via a potential V{r) that asymptotically acquires a 
dipole-dipole or van-der-Waals character, as sketched in 
Fig. [ija) . This distance dependence of the pair state po- 
tential curve is the main ingredient for generating a dis- 
tance selective loss feature. The energy of the pair state 
|ee) relative to \gg) is given by the two-atom detuning 
A2(r) = 2 A -I- V{r)/h. Thus, by adjusting A one can se- 
lect a critical distance at which the two-atom detuning 
vanishes (A2(rc) = 0) and Rydberg atom pairs are reso- 
nantly excited [T51[TH]. As sketched in Fig. [T|a), once a 
Rydberg atom pair is present two decay channels become 
available that ultimately lead to the loss of the excited 
atom pair. The atoms can either be expelled from the op- 
tical lattice by experiencing a repulsive mechanical force 
that accelerates the pair out of the system and/or by ac- 
quiring a strong momentum kick when spontaneously de- 
caying to a lower lying electronic state at (single-atom) 
decay rate F^yd. The combination of these two decay 
channels in conjunction with the distance-selective Ry- 
dberg excitation forms the basis of the non-local loss 
mechanism that eventually leads to dissipation-induced 
binding of atoms. 

With this qualitative picture in mind, let us now focus 
on the interplay between the non-local dissipation and 
the coherent hopping dynamics. The evolution of the 
system is governed by the Master equation p = — i [H, p\ + 
T^^p) with dissipator 
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The dissipative dynamics is characterized by the set of 
Lindblad jump operators Lj = y/j2a~ a~_^_^^, where 72 
corresponds to the effective decay rate of a pair of bosons 
populating two sites at the critical distance as depicted 
in Fig. [ijb). The parameter Ac = [^c/aj is an integer 
number measuring the critical distance in units of the 
lattice spacing a. 

Let us consider the regime of strong loss, i.e. 72 ^ J, 
in which we will now derive an effective equation. In this 
limit the dynamics is dominated by the dissipator ([iJ, 
whose steady state (eigenvalue zero) subspace is spanned 
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FIG. 2. Numerical propagation of two hard core bosons un- 
der the original master equation (with the dissipator . The 
graphs show the time evolution of the mean number of atoms 
on each lattice site. The initial state is cr^^crj |vac), the critical 
radius is Ac = 3 and the lattice has 20 sites and open bound- 
aries. Weak dissipation leads to quick loss of the particles 
from the lattice. Strong dissipation results in an overall re- 
duction of the decay and imposes a constraint on the coherent 
evolution that eflectuates a binding mechanism. 



by all arrangements of bosons in which there are no pairs 
occupying lattices sites at the distance Ac. The projec- 
tor on this subspace is given by Qo = Ylj (1 — ^jf^j+xj 
where rij = cr'^cfj . As sketched in Fig. [l|b) the co- 
herent Hamiltonian evolution can couple configurations 
from this subspace to the subspace containing one pair 
of bosons at the critical distance. Since 72 ^ J we can 
adiabatically eliminate these quickly decaying configura- 
tions [H] . This leads us to an effective Master equation for 
the projected density matrix p^s = QopQo with effective 
Hamiltonian 
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(2) 



and an effective dissipator that has the form ([!]) with the 
modified Lindblad operators 
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The effective Hamiltonian is simply the original hopping 
Hamiltonian projected on the subspace Qo- This projec- 
tion essentially amounts to a hard core interaction that 
forbids the simultaneous occupation of the sites j and 
j + Ac and ultimately is the origin of particle binding. 
Before discussing this in more detail, we briefly analyze 
the effective jump operators Lj^. An inspection of their 
structure shows that they lead to an incoherent decay of 
configurations that are "one hopping event away" from 
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a configuration that would decay under the fast dynam- 
ics with rate 72. An example of such a configuration 
is given in the upper panel of Fig. [l|b). More gener- 
ally, since the effective jump operators contain a sum of 
terms there exists in principle the possibility for certain 
superpositions of configurations of atoms not to decay 
although in each configuration the atoms are "one hop- 
ping event away" from being annihilated. When acting 
on these subradiant superpositions the contributions of 
the individual terms in the effective jump operators can- 
cel. However, this effect is not essential for the scenarios 
discussed in this work, as we focus on initial states that 
do not have subradiant contributions. The timescale of 
the effective dissipative dynamics is r oc 72 as can 
be read off from the prefactor of the effective jump oper- 
ators. Therefore in the limit of strong distance selective 
two-body loss (72 ^ J) dissipation is actually strongly 
suppressed and the dynamics of the system is dominated 
by the projected Hamiltonian HcS- 

To visualize this effect directly we study at first the 
evolution of a pair of bosons by numerically solving the 
original master equation for various loss rates. The criti- 
cal distance at which atoms decay is chosen to be Ac = 3 
and the initial conditions are such that two bosons oc- 
cupy adjacent lattice sites at the edge of the lattice, i.e., 
the initial state is cr^cr^ |vac), with |vac) being the empty 
lattice. Fig. [2] shows the evolution of the atomic density 
on a lattice with 20 sites and open boundaries, for four 
values of 72. For 72 = we find that after an initial pe- 
riod [Jt < 10) during which the wave packet stays local- 
ized and evolves ballistically, the density quickly becomes 
distributed over the entire lattice. For small decay, i.e. 
72 = J the short time behavior persists, but the decay 
leads to a quickly decreasing overall density. This situ- 
ation, however, changes if 72 is increased further. The 
wave packet remains localized, the propagation speed is 
approximately half of what is observed in the case 72 = 
and - most importantly - the overall density is much more 
slowly decaying as was observed in the case 72 = J. The 
latter is a direct consequence of the reduced decay rate 
connected the jump operators L'^^ that was discussed 
above. 

Let us now have a look at the spectrum of H^s in order 
to gain some analytical understanding of the binding as 
well as of the observed numerical results. We start with 
the case, in which the critical distance is Ac = 3 as in the 
above-discussed numerical example. Bound states occur 
once a number of particles is localized inside a region of 
the lattice whose extent is smaller than Ac. The distri- 
bution of particles cannot spread since doing so would 
imply the occupation of configurations that are annihi- 
lated under the action of the projector Qq. A possible 
dissipatively bound complex is, hence, formed by three 
bosons: \3j) ~ o-^^cr^(T+^^ |vac). This is a trivial so- 
lution of -ffoff with eigenenergy zero. A far more inter- 
esting scenario, which reveals that the emerging clusters 



may have an internal structure, is achieved for two par- 
ticles. In this case, analytic solutions of Hes can be ob- 
tained using the basis states |j,t) = (J^_i(j~^ |vac) and 
I:/, I) = CTjlj^iT^]^ |vac) which can be regarded as internal 
states of a two- atom complex localized at site j. The 
dynamics is then described by the effective Hamiltonian 
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which clearly displays spin-orbit coupling, i.e. the inter- 
nal "spin" state is coupled to the motional state. This 
becomes also evident from the eigenstates which are ob- 
tained via a Fourier transform: 
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As can be seen from this expression the internal state of 
the two-atom cluster (given in the brackets) indeed de- 
pends on the external quasi momentum qK = The 
dispersion relation of the complex is given by e{K)± — 
±2Jcos(2ii) which explains the difference in propaga- 
tion speed cx dx^iK) that was observed in Fig. [2} In the 
uppermost panel the bosons are approximately free (ne- 
glecting the hard core interaction) and the corresponding 
dispersion relation is oc cos {qk)- This leads to a maxi- 
mum propagation velocity for individual bosons that is 
twice as fast as that of the complex - matching the nu- 
merical data. 

Other situations with larger Ac and/or more parti- 
cles within the critical distance can be analyzed in a 
similar fashion. For Ac = 4 and three atoms localized 
within the critical distance a complex centered at site 
j has four internal states: | j, 1) = cr^^cr^cr^^ |vac), 
b',2) = o-+_2cr+o-+^i |vac), |j,3) = (T+_2cr+_iO-++i |vac) 
and |j, 4) = aj"_2<Jj'_^aj' \va,c) and its lattice dynamics 
is governed by the dispersion relation 



a + /3-\/9 -I- 8acos(qK) 



(5) 



which, opposed to the previous case, has four branches la- 
beled hy a, f3 — ±1. These few examples already demon- 
strate that strong non-local dissipation gives rise to a 
whole plethora of bound complexes with varying internal 
structure. 

Let us finally return in more detail to the question 
of how to realize the proposed scenario by exciting Ryd- 
berg atoms in optical lattices [21^ . In order to exclusively 
excite Rydberg pairs and not the single Rydberg atom 
states \eg) and \ge) the laser has to be far-detuned from 
the single atom transition, i.e., |A| Specifically, 
one needs to ensure that the photon scattering rate from 
these single Rydberg atom states 71 = r^ydl^^/A^ - with 
FRyd denoting the radiative decay rate of the Rydberg 
state - is much smaller than the hopping rate J. With 
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this being the case one may adiabatically ehminate the 
states \eg) and \ge) from the electronic dynamics. The re- 
maining pair states \gg) and |ee) form a two- level system 
with the Hamiltonian H^x — — ricff -I- \gg){ee\) + 

A2(r)|ee)(ee| where ilcft — r2^/(2A) is the effective Rabi 
frequency. The diagonal elements of this Hamiltonian 
represent the lowest and uppermost curve in Fig. [ija). 
Population is transferred between these potential curves 
when the wave packet of the relative motion is localized 
in the vicinity of the critical distance r^, as shown in Fig. 
[ija). Once transferred to the uppermost curve the atoms 
experience a strong mechanical force, F = —W{r) and 
are prone to radiative decay at (single-atom) rate Ffiyd. 
This excitation dynamics is completely analogous to that 
encountered in molecular photodissociation so that we 
can use the theory for diatomic molecules to estimate 
the atomic loss rate. On pair resonance (A2(rc) = 0) it 
is given by [2l] 72 = "i^l^f Re dt e~'^^^y''C{t), with 
the autocorrelation function C{t) = dr x* [i"-, 0)x('', t) 
and the wave packet of the relative motion x(r, t). In or- 
der to obtain an analytic expression for C{t) we employ 
the limit of a deep lattice potential and approximate the 
Wannier function at each lattice site by a Gaussian. The 
relative width a/a = [\/27r(Vo/i^r)^^'*]~^ of the Gaus- 
sian is controlled by the lattice potential depth Vo mea- 
sured in units of the recoil energy = n'^li} / {2ma?) 
(m denotes the mass of an atom) [22]. For Tc > a the 
overlap of the Wannier functions is zero, so that the 
wavepacket x(r, 0) describing the initial atomic pair at 
distance r > a is again a Gaussian with width ^/2a. To 
obtain a simple expression for C(t) we make a linear ap- 
proximation of the interaction potential around the criti- 
cal distance, V{r) « Vc—F{r—rc), with = drV{r)\r=r^- 
For short times we can neglect the kinetic energy of 
the atom pair and the wave packet evolves according to 
X(r,t) sa e-xjp[-iVct/h + iF{r - rc)t/h]x{r,0). Within 
these approximations the pair loss rate is 
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where erf(-) denotes the error function. Note that this 
rate is inversely proportional to the magnitude F of the 
force at the resonance point r^, i.e., the Rydberg pair 
excitation rate decreases with increasing mechanical force 
between the particles. The parameter regime, in which 
we expect non-local dissipation to form observable bound 
complexes is characterized by the inequalities 



71 < |J| <72- 



(7) 



In practice, the hopping rate J can be adjusted inde- 
pendently of 71 and 72 by changing the depth of the 
optical lattice potential. The single atom and two-atom 
excitation rates both depend on the Rabi frequency, the 
laser detuning and the radiative decay rate of the Ry- 
dberg state. However, since 71 is independent of the 



Rydberg-Rydberg interaction potential, it is possible to 
vary both rates independently by changing the force F 
at the resonance point. This potential shaping can, e.g., 
be done by applying weak external electric fields and 
thereby making use of the huge polarizability of Ryd- 
berg states [531 HI] • For example, using Rubidium atoms 
excited to the 455* Rydberg state with laser parame- 
ters (fJ, A) = (0.5, -30) MHz and adjusting F such that 
Fa/h — 2Tr X 100 kHz in an optical lattice of spacing 
a = 500 nm and depth Vq = lOE^ at a critical radius 
of Tc = 3a one obtains (71, J, 72) ~ (3, 50, 570) Hz, i.e., 
well within the regime in which dissipatively bound com- 
plexes are expected to form. With these parameters their 
lifetime is r w 0.2 s. 

In conclusion, we have studied the dynamics of a one- 
dimensional lattice gas of hard-core bosons with distance 
selective two-body loss. We have shown that this unusual 
dissipation leads to the formation of bound particle com- 
plexes, provided that the loss rate largely exceeds the 
coherent tunneling rate. Beyond a full classification of 
these dissipatively bound complexes in one dimension, it 
will be interesting to consider in the future higher di- 
mensional systems as one can expect here a much richer 
internal structure and more intricate dispersion relations. 
Moreover, since Rydberg atoms permit the realization of 
two-body interactions with angular dependence one can 
envisage a scenario in which dissipation does not depend 
only on the relative distance but also on the relative ori- 
entation of particles. 
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